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ABSTRACT 

According to the recycling scenario, millisecond pulsars (MSPs) have evolved 
from low-mass X-ray binaries (LMXBs). Their orbits are expected to be circular 
due to tidal interactions during the binary evolution, as observed in most of 
the binary MSPs. There are some peculiar systems that do not fit this picture. 

Three recent examples are PSRs J2234-I-06, J1946-I-3417 and J1950-I-2414, all 
of which are MSPs in eccentric orbits but with mass functions compatible with 
expected He white dwarf companions. It has been suggested these MSPs may 
have formed from delayed accretion-induced collapse of massive white dwarfs, 
or the eccentricity may be induced by dynamical interaction between the binary 
and a circumbinary disk. Assuming that the core density of accreting neutron 
stars in LMXBs may reach the density of quark deconhnement, which can lead 
to phase transition from neutron stars to strange quark stars, we show that the 
resultant MSPs are likely to have an eccentric orbit, due to the sudden loss of the 
gravitational mass of the neutron star during the transition. The eccentricities 
can be reproduced with a reasonable estimate of the mass loss. This scenario 
might also account for the formation of the youngest known X-ray binary Cir 
X—1, which also possesses a low-field compact star in an eccentric orbit. 

Subject headings: stars: neutron - X-rays: binaries - pulsars: individual (PSR 
J2234+06, PSR J1946+3417) 


2 


Introduction 


10^—10® G), rapidly rotating neutron 


stars (NSs) with spin p 

eriods < 30 ms ( 

Lorimer 

2008 

1 . In the standard recvcling see- 

nario ( 

Alpar et al. 

1982|: 

Radhakrishnan & Srinivasan 

1982 

Bhattacharva & van den Heuvel 


199lh . binary MSPs are regarded as the descendants of low-mass X-ray binaries (LMXBs), 
in which a NS has accreted sufficient mass from its companion star, and spun up to mil- 
l isecond periods. Rec ent observation s of the LMXB/MS P transition in PS Rs J1023-|-0038 
( Archibald et alj|2009 '). J1824—24521 f Papitto et ah 2013 ). and J1227—4853 (jRov et ahll2015 ) 
have lent strong support to the evolutionary link between LMXBs and MSPs. 

I f the initial orbital perio d of a LMXB is longer than the so-called bifurcation period 
(e.g., IPvlvser fc Savoniiel 119891) . the mass transfer begins when the companion star evolves 
to be a (sub)giant. The binary orbit continues expanding until the donor star loses most of 
its hydrogen envelope, forming a low-mass He white dwarf (WD). Since the mass transfer 
timescale is much l onger than tida l circularization timescale, the orbits of binary MSPs are 
highly circularized (IPhinnevi 119921) . except those in globular clusters where the orbits may 
be perturbed by dynamical interactions in the dense stellar environment. 

The discovery of PSR J1903-I-0327 presents a chall enge to the recycling sc enario, which 
is a binary MSP with a high eccentricity (e ~ 0.44) (jChampion et al.N2008f) . Located in 
the Galactic held, this pulsar is accompanied by a G-type main-sequence secondary. It was 


suggested to stem from a hierarchical triple system flFreire et al.ll201ll:lPortegies Zwart et ah 


2011 


Piiloo et al.j2012 ). or be a newborn NS experienced supernova fallback disk accre¬ 


tion ( Liu fc lJ 2009 ). Since then, some m ore eccentric MSPs have been d iscovered, includ- 


i ng PSRs J2234-I-06 ( Deneva et al.N2013 ). J1946-I-3417 ( Barr er al.l 12013 ) and J1950-I-2414 


(iKnispel et al.ll2015fl lH. All of these objects are Galactic held MSPs with a He WD (of mass 
around 0.24Mq for a 1.4Mq NS) in orbits with periods ranging from 22 to 32 days. The 
short spin periods ( 2 < Pg < 12 ms) indicate that they have experienced extensive mass 
transfer processes during the LMXB phase, but the eccentric orbits (with 0.027 < e < 0.13) 
are difficult to understand in standard LMXB evolution scenarios. Like in the case of PSR 
J1903-I-0327, it is still in principle possible that all these systems formed by the disruption 
of a triple system. However, the strong similarity in the orbital properties of all the recently 
discovered pulsars is not what one would expect from the disruption of a triple system: 
this is a chaotic process that should produce a wide range of orbital periods and orbital 


^Another MSP with an eccentric orbit will be presented in Camilo et al. (2015, in preparation). A 
previously known system, P SR J1618—3 919 might also be related: it has a spin period of 12 ms and an 
orbital eccentricity of 0.027 ( BailesI 2010 ). 
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eccentricities. The close similarity of these systems suggests instead a common evolutionary 
mechanism, different from what formed PSR J1903+0327. 


Freire &: TaurisI (120141) suggested that these MSPs formed from delayed accretion-induced 


collapse of massive WDs. According to their scenario, the progenitor binary may consist of 
a 7Mq primary and a 2Mq secondary. The primary star evolves to be an ONeMg WD of 
mass 1.2Mq while the secondary is still on its main sequence stage. In the following evo¬ 
lution, the secondary fills its Roche lobe and transfers hydrogen-rich material to the WD, 
so that the mass of WD increases to the Chandrasekhar limit. Howe ver, if the W D is fast 


spinning, it may keep accreting without collapsing to be a NS (e.g. lYoon fc Langerl 120041) . 


The WD-NS conversion only occurs when the WD has spun down and the centrifugal forces 
that sustain the star are weakened, which might take several Gyr. Then the sudden released 
gravitational binding energy during the collapse increases the orbital period and imposes an 
eccentricity to the former circularized orbit. The re-cir cularization times for such MSP/WD 
binaries are gener ally much longer th an Hubble time (IZahnI 119771: iHutIIl98lh . According to 
the calculation bv lTauris et al.l (120131) . there is an optimal range of orbital periods where the 
mass transfer is both dynamically stable (i.e., without entering common envelope evolution ) 


and rapid enough for stable thermonuclear burning of Hydrogen (e.g., lYaron et al.N2005l) . 


This range of orbital periods goes from 10 to 60 days, and the middle of this range is right 
where the aforementioned systems are. 


Alternatively, lAntoniadisI (120141) proposed that the eccentricities of the two MSPs may 
be caused by dynamical interaction between the binary and a circumbinary disk, which 
formed from the material escaping the donor star during hydrogen-shell flas hes shortly be¬ 


fore th e WD cooling phase. Adopting a linear perturbation theory following iDermine et ah 


(120131) . the author argued that, a disk with a fine-tuned lifetime as long as ~ 10^ yr and 
mass around a few IO^^Mq may result in eccentricities of e ~ 0.01 — 0.15 for orbital periods 
between 15 and 50 days. 


Adopting the idea that the orbital eccentricity is caused by sudden mass loss, here we 
suggest another possible scenario to this problem. During the recycling phase the NS accretes 
mass and angular momentum from the donor star. When its central density goes up to the 
critical density of quark deconfinement, a phase transition from a NS to a strange quark star 
(SS) may take place within a short time, and the loss of binding energy of the NS, if not big 
enough to disrupt the binary, can induce an eccentricity of the orbit. If the phase transition 
occurs at the end of or after the recycling process, the SS appears as an eccentric binary 
MSP with a He WD companion. 
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The NS-SS transition scenario 


Below we describe in some detail how the the NS-SS transition can influence the forma¬ 
tion of various types of MSPs. 


The concept of SSs has been suggested decades ago ( 

Itoh 

1970; 

Bodmer 

1971; 

Witten 

1984; 

Farhi & Jaffe 

1984: Haensel et ah 

1986; 

Alcock et ah 

1986 

), and it was suggested that. 


due to the absolute stability of quark matter, at least some of the pulsars could be SSs. 
Phase transition occurs when the central density of a NS rises above the critical density 
for dissolution of baryons into their quark constituents. The detailed processes of quark 
deconfl nement in massive NSs are not precisely kn own. It might be gradual, lasting around 
10® yr fjOlintol 119871: iHorvath fc Benvenutol Il988l) . or in a detonation mode, accompanied 


with energy releas ed compatible with a core collapse supernova (e.g., ICheng &: Pail Il996 


Ouved et ah 1200^ . Based on the standard equation of state (EOS) of neutron-rich matter. 
Staff et ahl (120061) pointed out that the critical density of quark deconflnement is 5po, 
where po ~ 2.7 x lO^^g cm“® is the nuclear saturation density. Although different EOSs give 
different masses when the critical density is reached, a gravi tational mass of l.S M^ may be 
a reasonable estimate for a NS with a core of quark matter (lAkmal et al.lll998h . 


Quite a f ew specific SS candidates (e.g., SAX J1808.8—3658. Ihi et al.lllOQOl) have been 


proposed (see IWeberl 120051. for a review), usually based on their observational features in 


X-rays. The formation of a SS can also impact the dynamical evolution of the host binary, 
which is seldom mentioned in the literature. We note that efficient mass growth of NSs is 
most likely to take place in LMXB^, and the phase transition may occur in the following 
two situations. 

In the first one, an accreting NS converts to be a SS during the LMXB stage when it 
becomes sufficiently massive. Some baryonic mass is abruptly transformed into its binding 
energy during the phase transition, and the orbit expands and becomes eccentric. The donor 
star is then detached from its Roche-lobe, temporarily halting the mass transfer, and the SS 
becomes a MSP. If the orbit is close enough, its high-energy radiation and particles might 


be able to ablate the companion star, leading to the formation of redbacks flRobertsI 12013 
for a review). Otherwise, after some time the companion star evolves to refill its Roche lobe 
and resumes the mass transfer. The orbit is then invariably circularized by tidal torques, so 
the binary evolves as a normal LMXB, producing a circular MSP/He WD binary. 


^NSs in intermediate-mass X-ray binaries usually accretes at a super-Eddington rate, so that most of 
the mass may be ejected. During the subsequent LMXB evolution when the accretion rate becomes sub- 
Eddington, the NSs can keep the accreted material and grow in mass. 
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In the second, if the NS’s spin is accelerated to be too fast, the central density might 
be centrifngally diluted and lie below the critical density for pure quark matter. Phase 
transition is delayed, and takes place during the terminal stages of the mass transfer or 
even the su bsequent MS P phase, during which the NS is spun down due to the propeller 
mechanism (iTaurisI 1201211 . or magnetic dipole and/or gravitational radiation, respectively. 
The binary will keep its eccentricity because of the extremely long circularization time for 
two compact component stars. That may be responsible for the the eccentric orbits of PSRs 
J2234+06 and J1946+3417. 


3. Phase transition-induced orbital change 


Similar as the lFreire fc TaurisI (1201411 model, here the orbital change in the MSP binaries 
is assumed to be due to instantaneous mass loss. To find out the possible parameter space 
of the WD mass and the mass loss during the phase transition for PSRs J2234+06 and 
J1946+3417 from the MSP mass, the orbital period, and the eccentricity, we use a Monte 
Carlo method to simulate the orbital change in the following steps. 

1. We assume that t he current mass Mss of t he MSP is randomly distributed in the 
range of 1.4 — 2 .OM 0 (IGangopadhyay et al.N2ni31 1. The upper limit is set in accord with 
recent measurements o f two massive pulsars, PSR J1614—2230 of ma ss 1.97(±0.04)Mfn 
( Dernorest et ah 2010 1 and PSR J0348+0432 of mass 2.01(±0.04)Mq ( Antoniadis et ah 

2OI3I 1. 


2. The current mass Mwd of the He WD is then calculated from Mss and the observ ed mass 
function M-^psin^f/(Mss + Mwd)^ = 0.0027 for PSR J2234+06 ( Deneva et ah 2013 1. where 
i is the inclination angle of the orbit. We assume that cosi is uniformly distributed between 
0 and 1. 

3. Assuming that the binary orbit is initially circular and the current eccentricity (assumed 
to be in the range of 0.11 — 0.15) is induced by instantaneous gravitational mass loss of the 
NS, we can get the amount of the ma ss loss AM = e • (Mwn + Mss) and th e gravitational 
mass Mns (= Mss + AM) of the NS (IBhattacharva fc van den Heuvellll99ll) . if the energy 
releasing during the phase transition is spherically symmetric and there is no kick imparted 
to the SS. 


4. Not all of the obtained values of AM and Mwd are physically reasonable, and we select 
those which satisfy the following condition predicted by LMXB evolution. On one hand, 
from Mss, Afwo and the current orbital period (assumed to be in the range of 22 — 32 days), 
we get the orbital separation a. Combining a and the eccentricity e yields the separation 
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before the phase transition via the equation Oj = a(l — e). On the other hand, we derive the 
orbital period Porb, i be fore the phase tran s ition from Mwd, using the orbital period - WD 
mass relation given bv iTauris &: Savoniiel (119991) for solar metallicities {Z = 0.02). Then 
from Porb, i, uud Mwd we also obtain Oj from the Kepler’s law. We require that the 

value of this Oi should match the former one within 1%. 

The simulated AM vs. M^d distribution is shown in Fig. 1. The allowed maximum 
NS mass is around 2.3Mq. It is seen that the mass loss is about 0.20 — O.34M0, and the WD 
mass is constrained to be about 0.257 — O.275M0 (red dots). If we decrease e down to 0.027 
(the eccentricity of PSR J1618—3919, which might have the same formation mechanism), 
the mass loss can be as low as about O.O5M0, while the upper limit of the possible WD mass 
increases to about 0.278(blue dots). Simulation with the mass function = 0.0039 for 
PSR J1946+3417 ( Barr er ah 20131) gives similar result, as shown in Fig. 2. These numbers 
are compatible with the standard evolution of LMXBs. 

The observed range of the orbital periods for eccentric MSPs thus far is 22 — 32 days, 
and the corresponding orbital periods before phase transition are found to range from 15 to 
35 days. It is not clear why phase transition would only occur for such a particular range of 
orbital periods. It might be related to the requirement that phase transition (supposedly in a 
massive NS) should occur at (or near) the en d of mass transfer, in or der to produce observable 
eccentric MSPs. According to the studies of lLi fc Wand fjl998l) and iTauris fc Savoniiel fjl999l) 
(see also Zhu et ah 2013), the amount of mass accreted by a NS during the LMXB evolution 
generally decreases wit h the orbital period, probably due to the thermal-viscous instability 
in the accretion disks fjLasotall200lh . In systems with longer orbital periods, the NS may 
not accrete sufficient mass to trigger phase transition, while NSs in narrower orbits may get 
enough material and undergo phase transition during the LMXB phase, but the resulting 
eccentric orbits would be circularized by the subsequent mass transfer. So we tentatively 
speculate that only LMXBs with a particular range of orbital periods may lead to the 
formation of eccentric MSPs. However, we need to caution that the accretion processes in 
LMXBs are rather complicated. Considering the fact the mass transfer is very likely to 
be nonconservative and the big uncertainties in the accretion efficiency, it is premature to 
accurately estimate the accreted mass by a NS. 


4. Discussion 

Based on the assumption that massive NSs may convert to be SSs, we suggest that the 
sudden loss of gravitational mass during the phase transition may account for the eccentric 
orbits in the two recycled binary pulsars, which are otherwise thought to be in circular orbits. 
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Recent population synthesis calculation bv IZhu et al.l (120131) suggested that about 0.1 — 
10% of LMXBs can produce SSs, depending on the masses of the nascent NSs and the fraction 
of transferred matter accreted by the NSs. Observational clues have also been proposed for 
the existence of SSs. For example, it is known that special global oscillation modes (r-modes) 
that emit gravitational waves would prohibit pulsars from spinning with high frequencies, 
unless the dampi ng of these modes, determined by the microscopic properti es of matter, can 
prevent this fe.g.. lAndersson fc Kokkotasll200lh . lAlford fc Schwenzerl (120141 1 showed that for 
each form of matter there is a distinct region in the spin frequency/spindown-rate diagram 
where r-modes can be present. They found that stars containing ungapped quark matter 
are consistent with both the observed radio and X-ray data, whereas, even when taking 
into account the considerable uncertainties, NS models with standard viscous damping are 
inconsistent with both data sets and additional damping mechanisms would be required. 

Our phase transition scenario can be compared with others in the following aspects, 
which are to be tested by future observations. 

1. Although phase transitions are expected to occur in massive NSs, th e masses of the 
resultant SS MSPs could lie in a wide range (IGangopadhvav et ahl l2013l) . depending on 
the amount of the m ass that is transferred into the binding energy. This is similar to 
the Antonia^ ( 2014 ') prediction but different from the simplest 


scenario m 


Freire fc Tauris 


(120141) (that of a rigidly rotating WD before the collapse), which predicts that these MSPs 
should be exceptionally light. 

2. The eccentric MSPs probably have a He WD companion with mass in a relatively narrow 
range (~ 0.25 — 0.28 Mq). With such WD companions, the NSs are likely to be heavy at 
the moment of phase transition, implying that the resultant SS pulsars are fully recycled 
pulsars. 

3. The eccentric MSPs may have kinematic properties (i.e.. Galactic heights and peculiar 
velocities) at least similar as what is observed in the normal MSP population, since the 
super nova that formed t he NS should give the system a reason able kick. This i s again similar 
to the lAntoniadisI (12014J ) prediction but contrasts with that of iFreire fc TanrisI (120141 ). In the 
latter scenario, the AIG is the only supernova that ever occurs and, because this supernova 
is so symmetric, the eccentric systems generated would have very small peculiar velocities 
and Galactic heights. 

As noted before, phase transition may also occ ur during the LMXB pha se. The con¬ 
version of a NS into a SS could lead to a quark nova (jOuved et ahl l2002l . 120111) . with energy 
released comparable to a normal supernova. A peculiar X-ray bin ary, Gir X—1, might be p ro- 
duced in this way. This source has an orbital period of 16.6 days (iKaluzienski et al.lll976l) . It 







































does no t seem to fit either a high-mass or a low-m ass X-ray binary. Detection of type-I X-ray 
bursts flTennant et al.l Il986l: iLinares et al.l l2010h firmly establishes its nature a.s a low-field 
(< 1 0 ^° G) NS, s i milar as those in LMXBs. The observed spectr al features flShirey et al. 


1999 : Soleri et al. 2009 ) and kilo-Hertz quasi-periodic oscillations ( Boutloukos et al. 20061) 


also lend strong support to t his conjecture. However, the orbit is shown to be eccentric with 
e = 0.45 (jjonker et al.l 120071) . suggesting that the binary was recently s ubject to mass loss 
(and/or a kick) during the formation of the compact star. Moreover, iHeinz et abl (120131) 
discovered a supernova remnant (with age less than 4600 y r) surrounding Cir X—1, making 


it the youngest known X-ray binary. Honker et al.l (120071) associated the observed optical 


absorption-line spectrum with the companion star of Cir X—1, suggesting that the compan¬ 
ion star could be a B5—AO supergiant. However, the possibility cannot be excluded that 
the spectrum originates in the accretion disk, and in that scenario, the companion star is 
likely to be a (sub)giant of mass around O.4M0. If it is the case the conflicting properties 
(young age, eccentric orbit, and low magnetic held) of Cir X—1 can be well accounted for 
by the NS-SS phase transition scenario. The key point is that the supernova remnant and 
the eccentric orbits were produced by the quark nova when the accreting NS was converted 
to be a SS, rather than the supernova that originally formed the NS. 


No matter if our specihc scenario is conhrmed or ruled out by future observations, 
the present data discussed in this paper will provide useful constraints on the formation 
channels of these peculiar systems, and there will be deep consequences for the physics of 
strong interactions. 
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White dwarf mass (Mq) 


Fig. 1.— Required gravitational mass loss of the NS vs. the WD mass in the phase transition 
scenario to reproduce the properties of PSR J2234+06. The red and bine dots represent the 
cases for the eccentricity in the range of 0.11 — 0.15 and 0.027 — 0.11, respectively. 
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White dwarf mass (Mq) 


Fig. 2.— Same as Fig. 1 but with the mass function for PSR J1946+3417 



